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ABSTRACT 


The  problem  of  predicting  the  induced  roll  phe¬ 
nomena  associated  with  finned  missiles  at  high  angles  of 
attack  has  been  approached  by  examining  the  interaction  be¬ 
tween  the  body  wake  vortex  system  and  the  finned  surfaces 
of  the  missile.  The  application  of  the  classic  Blasius 
Moment  Integral  to  this  problem  permits  the  independent 
evaluation  of  each  singularity  in  the  flow  field  which  may 
contribute  to  the  moments.  This  method  is  applicable  to 
missiles  employing  fin  schemes  of  any  geometry  and  is  used 
herein  to  analyze  the  induced  roll  of  a  typical  cruciform 
missile  configuration. 

A  computer  program  has  been  developed  to  generate 
numerical  solutions  utilizing  those  parameters  which  affect 
the  magnitude  and  direction  of  the  induced  roll. 
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SECTION  I. 


Introduction 


Rockets  and  missiles  employing  fins  of  cruciform 
arrangement  have  teen  found  to  develop  anomalous  aerodynam¬ 
ic  rolling  moments  which  are  functions  of,  among  other 
things,  the  angle  of  attack,  the  bank  orientation  of  the 
wings  and  Mach  number.  These  rolling  moments  tend  to  up¬ 
set  the  stability  by  rolling  the  missile,  and  unless  they 
are  counteracted  by  roll-control  measures,  the  actual  flight 
path  control  accuracy  is  apt  to  deteriorate  at  a  crucial 
time.  It  is  generally  accepted  that  the  absolute  magnitude 
of  the  induced  rolling  moment  is  of  primary  concern. 

Induced  roll  phenomena,  which  could  be  significant 
for  some  configurations  at  low  angles  of  attack,  become  of 
increased  importance  as  the  angle-of-attack  range  is  extend¬ 
ed  for  the  purpose  of  increased  maneuverability.  The  source 
of  high  angle  of  attack  problems  has  been  found  to  be  as¬ 
sociated  with  body-wing  or  body-tail  interference. 

Without  the  mitigation  of  these  induced  rolling 
moments,  high-angle-of-attack  maneuvers  in  some  cases  may 
be  precluded  or  in  other  cases  result  in  increased  control 
system  complexity  or  increased  size  of  the  wing  surfaces. 
Such  alternatives  usually  result  in  overall  performance 
degradations. 

The  possible  major  sources  of  anomalous  induced 
rolling  moments  for  typical  rocket  configurations  were 
discussed  and  evaluated  in  Ref.  1.  The  magnitude  of  some 
of  these  rolling  moments  have  been  found  to  be  negligible 
compared  to  the  induced  rolling  moments  observed  from  ex¬ 
perimental  investigations.  Strong  evidence  is  available 
(Ref.  2)  to  support  the  argument  that  the  interaction  of 
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the  wing  or  tail  surfaces  with  the  vortex  system  generated 
by  the  body  is  largely  responsible  for  the  induced  rolling 
moment  phenomenon. 

The  vortices  formed  in  the  wake  of  a  body  affect 
the  attached  lifting  surfaces  by  altering  the  local  flow 
characteristics — i.e.,  angle  of  attack,  Mach  number,  static 
pressure,  etc. — across  the  span.  The  variations  in  the 
magnitude  and  spanwise  distribution  of  the  panel  normal 
forces,  which  are  produced  by  these  changes  in  local  flow 
characteristics,  are,  therefore,  functions  of  both  the 
strength  of  the  vortices  and  the  geometry  of  the  wing- 
vortex  system.  Consequently,  comprehension  and  prediction 
of  the  effects  of  wing-vortex  interactions  must  depend 
critically  on  the  description  of  the  vortex  system  itself. 
Fortunately,  experimental  data  are  available  (Ref.  3'  which 
describe  the  wake  vortex  system  generated  by  tile  body  of 
the  wing -body  configuration  considered  in  the  present  study. 
These  vortex  data  eliminate  the  additional  complexity  of 
relying  on  theoretical  estimates  of  the  vortex  system  and 
are  used  herein  as  a  basis  for  analysis. 

It  should  be  noted  that  the  Blasius  moment  inte¬ 
gral  has  been  recently  applied  (see  Ref.  I)  in  evaluating 
the  induced  roll  of  planar  missile  configurations. 
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SECTION  II.  General  Formulation  of  the  Problem 

In  order  to  simplify  the  problem  we  will  assume 
(as  was  done  in  Ref a  1  for  the  planar  model)  that  the  flow 
is  steady;  incompressible,  and  inviscid.  The  assumption  of 
negligible  missile  roll  rate  will  allow  for  the  steady  flow 
condition . 

A  typical  body-tail  configuration  utilizing  four 
equally-spaced  rectangular  fins  will  serve  as  our  model  for 
analysis;  (as  will  be  shown  later  the  fin  can  be  of  any  de¬ 
sired  shape  such  as  triangular,  trapezoidal,  etc.).  The 
model  in  the  cross-flow  plane  (Fig.  I)  shows  the  cross  flow 
velocity,  XJ*  /  the  missile  roll  attitude,  <f>  ,  and  pertinent 
model  dimensions.  IV/ 


VooSlV  Of — XT 


Figure  1.  The  Cross-Flow  Plane 

We  may  now  represent  the  complex  potential  for 
this  flow  in  the  Z  plane  (Fig.  2)  by  W”(Z  )  where  the  two 
body  vortices  are  shown  relative  to  the  fins  and  the 


cross-flow  velocity  TJ.  Because  of  the  missile  roll  the 
vortices  are  not  necessarily  symmetric  or  of  equal  strengths. 


\iy 


Figure  2.  The  Z  Plane 

For  steady  state  rolling  moments  about  the  origin 
the  Blasius  moment  integral  is  given  as 

where  the  contour  is  along  the  outside  of  the  body  surface. 

The  analysis  can  be  considerably  simplified  if 
the  configuration  in  the  Z  plane  is  conformally  mapped  into 
a  circle-  The  mapping  procedure  for  the  cruciform  shape 
can  be  found  in  Ref.  4.  Thus  the  plane  (Fig.  3)  repre¬ 
sents  the  mapping  of  the  missile  contour  into  a  circle  of 
radius  c .  (Note  that  the  magnitude  and  direction  of  the 


cross-flow  velocity  and  the  magnitude  of  the  two  body 
vortices  are  unaltered  by  the  mapping.*) 


Figure  3.  The  ^  plane 

The  relationships  between  corresponding  points  in  both 
planes  are 


z=± 


&*+  4s ) 


where  C 


-■\/2{bz+gn  A 


See  Appendix  A  of  Reference  1 


The  choice  of  signs  in  (2a)  and  (2b)  is  based  on  the  fol¬ 
lowing  rules: 


(JL)  The  "outside"  plus  and  minus  signs  hold  for 
the  right  half  and  the  left  half  of  the  Z 
or^  planes /  respectively. 

(2)  The  "inside"  plus  and  minus  signs  are  chosen 
accordingly  to  insure  that  the  exterior  of 
the  body  maps  to  the  exterior  of  the  circle. 


The  reader  is  referred  to  Ref.  1  (p.  15)  for  the 
derivation  of  the  complex  potential  w(^)  which  was  ob¬ 
tained  in  the  following  form: 


W(¥)  =  U  (elcp\  +  eLcp^) 

~  iVz  1 


4-  ifi 
2TT 


ZVr 


ZTT  C 


where  the  velocity  potential  takes  into  account  the  pos¬ 
sible  presence  of  a  Kutta  condition  by  including  an  arbi¬ 
trary  votex  of  strength  with  center  at  the  origin. 

As  shown  in  Ref.  1  (p.  16)  the  Blasius  integral 
may  also  be  expressed  as 


m'=  j>  /C?)  w Y?)  d  § 


Z 


where 
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and  W(fj  H 


The  use  of  the  expression  for  2  as  a  function  of  ^  (Equ.  (2b)} 
and  its  first  derivative 

dz  i  t-„  c*\  ■y(f+cVs*')>--‘iKL> 

d?  V?v  &)  e* 


gj.ves,  after  algebraic  expansion, 
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C^+c^-cD 


s%c4  + 

"VT?‘t+C‘*)-^4S4  -<'5*+c^ 


The  reader  is.  referred  to  Appendix  A  for  a  list  of  deriva¬ 
tives  of  ► 

The  expression  for  W2ft£)  can  be  found  in  Ref.  1 
(p.  17)  where  the  singularities  are  shown  to  be  poles, 
located  at 


1  =  0,'?,,^, 
c?  cz 


(The  poles  at  ^  and  are  inside  the  circle.)  Since 
and  *5*^ are-  exterior  to  the  circle^  ~C,  they  need  not 
be  accounted  for  as  contributions  to  the  integral. 

The-  remaining  singularities  of  the  integrand  all 
lie  on  the  circle  and  are  due  to^(^)  (Equ.  (5)).  These  are 
poles  at*'?-’  =  ±  c  and  ^  =  t(  c  and  branch  points  where  the 


(5) 


s 


Thus-  branch  points  are  located  at 

?  =  ±V£ 

where 


That  is,  eight  branch  points  appear  at 


(6) 
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These  are  all  on  the  circle  since  in  each  case 

\ 

|§j=  (C^+a^f2  =  (c*)'h=  c  . 
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SECTION  III.  Contour  Intearal  Evaluation 


The  procedure  for  evaluation  of  the  poles  is 
shown  in  Ref.  1  (p.  20).  Thus,  if  a  complex  function  has 
a  pole  of  order  ir.  at^^  ,  then  its  residue  there  is 

P  =  (vbyi  f(§)  •  ( 

Referring  to  Ecu.  (4) ,  a  contour  around  a  pole  of  order  m 
of  the  Blasius  integral  located  at  can  be  expressed  as 

Xr  =  -J tetri  fpQs~,rri)  ( 

where  can  be  defined  as 

F(i)=j(%)wzC s)  ( 

in  Equ.  C7) . 

For  the  simple  poles  at  =c  the  integral  is 

„\  '5  W'W)  V>e*5v  \ 


-fw; 

c  /  £-*C  ‘-J  -C2X 


;3+C2)  + 


= -^e7Tt-^  Jc2-t  qz+cx  J  Wz(c)  . 

Similarly,  for  the  pole  at  ^  =-c  the  integral  is 

I-c=  -^4  W*a* ±c~]  • 


The  total  contribution  due  to  the  poles  at  c  and  -c  is  thus 
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a*±cz 


wKc)i-wH-c) 


(10) 


Since  Re  w(c)=ReVV£c)  ~  O  J  it  follows  that 
Rei^C^ReiW2^^  and  therefore 

l'c+lLc.=  0  . 

Similarly/  the  total  contribution  due  to  the  poles 
at.  Lc.  and  ~£c  is 


(11) 


^LC.^-LcT  feTTL  =£ 


C  + 


a*+c* 


w*(ic)  +  \t/*(-Lc) 


(12) 


Since  Inr  w(£c)  =  Ifn  w(-ic)=  O,  it  follows  that 
Re  tw2U'c)  =•  Re  £w^(-£c)  =  O  and  therefore 

l'ic+-I-ic  =  o  . 


(13) 


The  contribution  to  the  total  moment  due  to  the 
simple  pole  at  ~  O  is 


r:=-^fiw  [-  (|f  ( A 

Since  ^(0)  =  O  (Equ.  (5))/ 


(14) 

The  contribution  of  the  second  order  pole  at 
=  0  is  of  the  identical  form  as  Equ.  (21)  in  Ref.  1  where 
-/(0)=0  and  (0)=-l  (see  Appendix  A).  Thus,  after  differen- 
tation,  elimination  of  imaginary  numbers,  and  collection  of 
terms  we  have,  as  in  Equ.  (22)  in  Ref.  1, 


I'-Rfefue*  |n(f 


The  contributions  to  the  total  moment  due  to  the 
third  and  fourth  order  poles  at  tjf  =0  are 


1 3  =  o 

_L.0 


since  ^  (01=0  {.see  Appendix  A) ,  and 


r^=o 


since  (0)=0  {see  Appendix  A). 

For  the  simple  poles  at  §  =  and  §  =  the 
contribution  to  the  total  moment  is  of  the  identical  fo 


contribution  to  the  total  moment  is  of  the  identical  form 
as  Equs.  (251  and  (26)  of  Ref.  1,  respectively. 

_  Finally,  the  remaining  second  order  poles  at 
|f=  and  -=r  are  similarly  evaluated  as  Equ.  (27)  in 
Ref.  1.  Expressions  for  -ft  in  Equs.  (25)  and  (26)  and  for 
in  Equ.  (27)  are  contained  in  Appendix  A. 

The  remaining  singularities  to  be  evaluated  are 

the  eight  branch  points  * - '^>VIII^  wll*ck 

on  the  circle.  We  may  conveniently  choose  branch  cuts 
joining  these  branch  points  as  straight  line  segments  as 
shown  in  Fig.  4. 


Figure  4.  Branch  points  and  branch  lines  of  /(■§> 


Let  us  define  the  coordinates  J  and  K  as 


We  can  enclose  each  branch  line,  e.g.,  ^  x  to  Ij;  IX  et 
with  the  contour  ABCDEFA  as  shown  in  Figure  5. 


Figure  5.  Branch  Line  Contour 
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Referring  to  equations  (4)  and  (5) ,  we  wish  to 
evaluate  the  integral 


4BCDEFA 


The  evaluation  of  the  integrals  around  AEC  and 
DEF  was  shown  to  be  equal  to  zero  (see  Appendix  C  in  Ref.  1) . 

Using  the  analysis  contained  in  Ref.  1  for  the 
evaluation  of  the  integrals  along  the  lines  FA  and  DC,  we 
find  that 


On  FA, 

"5  *^zr  anc^  = 
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r  r2< 

I  =  +  - — (lg--Osw^i c-r")  jr 

Jr-.  '  r)s4- E=  ~ k'V 


[ftxc-^+cU^-ry-c^Vfzkr-^r.  (23) 


(2 .  x-  r)(2  ?3r-  r) — -  -  r+(J+ ;  K)-)  +  (( r)  t-C4)] 


Similarly,  by  extending  the  analysis  along  CD  we  find  that 


=  TRcZfK 


vf _ i 

J  KW-r 


tV2fetr~r) cjt 


k-rr+c' 


r-fpC*lf-rW2^-r)r* 


(2;J--r)(2?k-r) - fer-r+(J+i< 


where  the  only  change  in  the  integral  is  indicated  by 

Equ.  (21),  Addition  of  (23)  and  (24)  yields,  after  algebraic 

reduction 

2K 

T  =  f ,  f  +<&*£.**¥  f ,Tgh(r-)^^e)r 

J  j  '  J  (/[*-c*)[j?(ir) + 0^+c¥)2] 

FA  CO  O  u 


where 


f  (r)  s  tar-r  =(K-r)+-iJ 


1? 


and 

/l  =  jrr\j(7H-r)r(2i  J-f)(2K+2iX-  r )(K+-  3+ £  (<T-K)-r )  • 

121 1 


The  same  procedure  is  used  for  the  evaluation  of 
the-  integral  around  the  branch  cut  connection  ^  with 
and  produces 

•2K 

~r  f  h(r)W  (rl')dr _ 

^ar^ET  _  >  I  (^-C^h^+C^+C^)2] 


where 


»j(0  ~  ?nr""lr  “  (V“ -r )—i>T 


r)f2/c--2£Cr-r)(^  J- 


It  should  be  noted  that  the  only  difference  between  the 
integrals  of  (25)  and  (28)  is  in  the  sign  of  L J. 
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Fiorure  6.  Branch  Line  Contour 


The  integral  v/e  wish  to  evaluate  is 


-  -  -  -£-W)  i 

a  'a 


ascdef 


Now,  we  may  use  the  polar  coordinates 

t-WsrEete“ 

?-^3ZH=  Tjzne^331 


(36a) 


7T 


where  both  ©•  _  and  ©  _T  are  allowed  to  vary  from  —  ~  to 

3TT  V  * 

2 


(35) 


(36b) 
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On  and©VII=^- 

On  C  ‘  D  ‘ ,  ©  VI=  —■  and  OVII=  JL 


Thus, 

®32i  ®izrr  =  { 

Therefore , 


2 

q  on  F'A' 
^TT  on  C'D'  . 


VS-fsiYS-'ini)  =V^  e‘^ 


Qaa-f-Qvir 


CO 


= 


“\/  fe  ^3231  G  on  F*A 

-Y^ 


on  C  P 


From  (36b),  along  F’A*  and  C'D1 

^Ysar+isne1®™ 

=  '?-sm+  ftm  eilr/z 

—  ^3ZE  +  l  f3ZZI 
and  d^=  idr^r  , 

Also,  on  F'A1  we  note  that 

-  (im  +  i rim) 

=  2(j-LK')+ir 

where  r=rVII  and  similarly 


.  (37) 
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"*  J~ If'" —  2ZT ■hit' 
i  -  J^iK+ir-^kf+il)-  J  +K+  tr~  (.  (*T +K) 
J~ik'^L)r"(K+LT^~J--l<!+Lf--L  (j~ +%') 

%~^=T-ik:+Lr^K-LD)= J+K+ur+i(j’-ic)  B 


Also 


V r 

Substitution  of  (37)  into  (35)  gives,  on  F 5 A 1 

2K 

I  .  rp.  (tznL-h'''rfw'Z(^Tm+ir ^(idr) 


r)r- 


(38) 


(j~+K-t-ir-l(j+l<)) - (2T^lrj~  ±((§Mr+^')‘f+c‘fL)]  _ 


Similarly,  on  CD1,  where  rVII  varies  between  2K  and  0 

.O 

f  i! tY  r^gg-Kt^^t-d^ _ 

j  =  -/Ce  f  J 


(39) 


(j+K-t-ir-i(J+K)) - — (zJ +;r)  ± ((%m+ir)U c1')] 
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Adding  (38)  to  (39)  we  obtain 


f  f  T  _  -  p,  g,.  nc#  I  1 sP(r)W*(r)  idr 

J+  = 

^  rW  o 


F'A*  W 


where 


T  =  ^~vrr  ~F  ir  —  3~iK+  if . 


Since  n  =  +  L 


(see  Equ.  (26)) 


then 


and  therefore. 


Also, 


p(r)  s  ^J(2K-r)r(^^  i  (J-+lty(J^&ir-i  (T+fOV 


(j -K+ir+i(^ 

We  may  also  show  that  _ _ _ _ » 

ip  =  V(TF  p=-f^/(2^-r')r(i.('JVK')-r +(J‘+k))  - 

(iC^-^^Xif^-r+K-r;  • 


(see  Equ.  (27) ) . 
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Thus,  .  2 

+h =-p . 


C43) 


The  same  procedure  is  used  for  the  branch  cut 

connecting  v  with  ^  vill*  Replacing  T  by  — land  P  by 

— p  we  thus  obtain 
r  2K 

f  (+Ts)(+p)\*/2(~T)cJr 

_  .  jv  i.x-  .  - 1  y  - - - - - - 

o 


*  " +&‘hfZ  Jic+t2>c23[(+T2VC^tP^+ft4+c^' 


(44) 


If  we  let 


p„  ±rfp _ 

x  (T*+C2Xr2-Ca)[^p*+(T4+C4)2] 


then 

£  = 


+ife 


<t2+cj)(f2-c2)[-pa+(r^+c4X' 


Therefore 


2K 


[R^W+RwY-T)]dr 


(45) 


(46) 


(47) 


or,  alternatively 


?K 


ww  -  *  «f  *7 fK-fl-'W]* 


(48) 


We  have  proved  that 


T-  —in 

-pa=h* 


(see  Equ.  (41) ) 
(see  Equ .  (43 ) ) . 


and 
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Substitution  into  Equ.  (45)  gives 

P=  -ifVi _ 

r  f  r|  z+cOfif-C*^ ha+(>J,4  C*)5] 

-n5h  _  =_M 

(T *Ic*J0f+  ca)[h%(T * +  c*)2] 

and  therefore, 

ZK 

±fe<Kf' 

o 

Adding  this  result  to  Equ.f34)will  give  the  total  moment 
contributed  by  the  four  branch  lines 

'■^nr/sasr + “^si/^sff +  ~ 

ZK 

O 

-[w2(^>wa(V0]j<]r  • 

In  summary,  the  total  rolling  moment  is  the  sum  of  the 
contour  integrals  around  the  singularities 

M'o.**  +^°  +^° 

+■  Xc«/f(  +  ifcy ta + i^/f,  +■  JVfz 

'^X/'?2T^  ,^3Sr*'^~^hZL  • 


(49) 


(50) 


(51) 


(52) 


i 


> 
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The  equation  for  each  integral  in  (52)  is  conveniently 
tabulated  below. 


Table  1 
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The  final  result  presented  by  Equ.  (52)  is  pro¬ 
grammed  (see  Section  IV)  to  obtain  numerical  solutions  with 
the  use  of  a  high  speed  digital  computer.  It  is  apparent, 
however,  that  results  may  be  obtained  for  limiting  cases. 
Thus,  for  example,  if  we  consider  that  no  external  vortices 
are  assumed  to  be  present  in  the  flow  fields  ( fl =  f**=0)  , 
then  Equation  (15)  equals  zero  as  well  as  Equations  (25) , 
(26),  and  (27),  of  Ref.  1.  Therefore  the  total  sectional 
rolling  moment  in  (52)  is  reduced  to 

m'o  =  ±  H 


-  [wY q)-w 


Using  the  expression  for  in  Ref.  1  we  may  note  that 

wX-lrf)=uie^ 


uy 

7 T 


e^0  + 

n  n3 


and 

w%in)  =  u*yi<p+  -&+?*-£)  + 


irr 

7 r 


L  n 


Gi<p  ,yl<P  c* 

“  ^  »3 


7 


therefore. 


w =  U2  [(e*0-  eai0)+  ^(esi0-e*0)  ] 


(53) 


(54) 


=  -7U*L$ir)Z(p 


c* 


+ 


n  L  «f.  * 
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V 


By  the  same  procedure, 

W*(n)-V*(rd=  -*u*i sm 2<p[~ ^J+  [l-  f*  ]  •  (55) 


After  subtracting  (55} from  (54)  and  substituting  the  result 
in  (5  3)  we  obtain 

o 

<ZK 

=  x&VfB+TR  [wf  ('- y^ysm  <f>-  cos  <?>)]dr 


(56) 


The  expression  for  H  (obtained  from  Equ.  (49))  is  now  sub¬ 
stituted  in  (56)  to  obtain  the  final  result; 


r*K  -  ~ 


(57) 


(See  Appendix  B  for  a  derivation  of 

Equation  (57)  may  now  be  investigated  in  the  fol¬ 
lowing  for  two  dimensional  flow,  in  the  absence  of  external 
vortices,  past  certain  configurations  which  are  oriented 
at  a  roll  angle  0  with  respect  to  the  cross  flow  velocity 
vector  t J  . 

A  simple  cross  configuration  is  obtained  if  we 
reduce  the  body  radius  R  to  zero.  If  R=0,  then  it  is 
obvious  that  (5/)  equals  zero  and  that  this  configuration 
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experiences  zero  net  rolling  moment  at  any  roll  angle  , 
free  stream  flow  density  and  velocity  TJ  .*  yhe  other 

limiting  case  to  consider  is  that  of  a  finless  cylinder 
where  R=b  (see  Fig.  1) .  Because  c=b  (by  the  conformal  trans¬ 
formation  of  the  Z  plane  to  the  ^  plans)  it  is  obvious  that, 
since  K=0  (see  Equ.  (19)),  Equ.  (57)  equals  zero  since  the 
upper  limit  of  integration  vanishes.  Thus  the  rolling  moment 
M'0,  for  a  finless  cylinder  is  zero  for  this  case,  as  ex¬ 
pected. 


* 

See  Reference  5. 
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SECTION  IV.  Numerical  Analysis  and  Digital  Computer  Results 

The  sectional  rolling  moment  M0  has  been  found 
(see  Equo  (52)}  to  be  a  function  of  the  following  para¬ 
meters  : 

mL = m i  (b, ,igs ,  n ,  p2 , f , tr, 4>, r)  . 

As  is  shown  in  Ref.  1  (Equ.  (43)  )#  /Vf^may  -be  expressed  in 
dimensionless  form  where 

CmI  as Cm'0  )  -j=r  )  -j=r  )Cr,  )Cr2')(P)  zrrfzd^)  ' 

Equation  (58)  is  used  to  obtain  numerical  data  with  the  aid 
of  a  digital  computer.  For  presentation  of  the  computed 
data,  however,  the  following  changes  were  made  in  the  above 
parameters: 

(1)  The  reciprocal  form  b/R,  is  plotted  with  b 
used  as  the  independent  variable. 

(2)  Vortex  positions  and  are  presented 

as  N/z  and  fel/R  ,  respectively,  which 
are  obtained  by  mapping  the  positions  in  the 

plane  to  the  corresponding  positions  in 
the  'Z  plane. 

(3)  The  free  vortex  circulation,  Y ,  is  dependent 
on  the  presence  of  a  Kutta  condition  and  is 
therefore  not  specified  independently. 

The  two  dimensional  (rectangular)  fin  and  vortex 
data  used  for  the  wind  tunnel  model  in  Ref.  3  provided  the 
inputs  to  the  computer  program  to  generate  numerical  data. 
For  a  two  dimensional  fin,  the  total  moment  in  dimensionless 


(58) 
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form  is  defined  as 


__  AAo 
C-Mo~  -LfTJJb2L 

where  M  =  the  total  moment,  L  =  the  chord  length. 

The  wind  tunnel  model  geometry  and  flow  conditions  used  for 
the  parametric  study  are  as  follows: 
b=  2.8  in. 

L=  2.5  in. 

R=  .75  in. 

0=  430  in. /sec 

^=1.45  X  10  ^  slugs/in. ^ 

<X-  20  degrees 
F>fV  285  in.^/sec. 

-  1.312  in.  (Separated  by  an  included  angle  of 
36  degrees) 

(It  was.  reported  by  C.  Wong  (Ref.  3)  that  these  vortex 
positions,  are-  essentially  symmetric  with  respect  to  the 
body . ) 

To  illustrate  the  vortex-fin  interaction  at  various 
roll  angles  of  interest,  Fig.  7  shows  the  relative  positions 
of  the  vortices  with  respect  to  the  fins.  The  numerical 
results  obtained  from  the  computer  study  are  shown  in  Figures 
8-13. 

(To  check  the  validity  of  the  above  results  this 
writer  performed  a  hand  calculation  of  the  individual  rerms 
of  Equ.  (52).  Except  for  Equ.  (51),  which  proved  difficult 
and  time  consuming  to  hand  calculate,  the  results  agreed  to 
within  five  percent.) 
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SECTION  V.  Conclusions 


Induced  rolling  moments  at  high  angles  of  attack 
are  due  principally  to  the  interaction  of  the  vortex  pair 
generated  by  the  body  with  the  attached  fins.  It  is  evi¬ 
dent  that  these  rolling  moments  are  strong  nonlinear  functions 
of  fin  span  and  vortex  strengths  and  positions. 

One  important  assumption  made  for  this  report  is 
that  the  flow  is  two  dimensional.  If  this  restriction  is 
removed,  the  problem  is  further  complicated  in  that  wing 
tip  vortices  as  well  as  the  body  vortices  must  be  accounted 
for.  This  is  necessary  as  the  body  vortices  are  undoubted¬ 
ly  influenced  by  the  wing  tip  vortex  system  as  they  pass 
over  the  fins.  The  resulting  interaction  between  the  vor¬ 
tices  of  both  systems  depends  on  their  relative  positions 
in  the  flow,  relative  strengths  and  directions  of  rotation. 

For  example,  two  vortices  of  equal  rotational  directions 
will  induce  velocities  on  one  another  causing  them  to  move 
in  opposite  directions  ,  while  two  vortices  of  opposite  ro¬ 
tations  will  move  towards  each  other.  The  above  inter¬ 
actions,  combined  with  vortex-fin  interaction  and  missile 
roll  motions  disturb  the  flow  symmetry  which  would  other¬ 
wise  prevail  if  these  effects  did  not  take  place. 

The  vortex  data  used  for  this  report  was  obtained 
by  experiment  with  the  use  of  the  subsonic  wind  tunnel  at 
Boston  University's  College  of  Engineering.  The  preliminary 
data  (see  Section  IV)  indicates  that  the  body  vortices  are 
symmetric  with  respect  to  the  body,  even  in  the  vicinity 
of  the  model  fins.  This  observation  was  made  after  reviewing 
photographs  of  the  missile  model  and  vortex  geometry  n  the 
wind  tunnel. 
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The  variation  of  roll  moment  with  roll  angle  is 
shown  graphically  in  Figure  8.  It  is  observed  that,  as 
the  roll  angle  approaches  18  degrees,  vortex  2  approaches 
the  horizontal  fin  causing  a  rapid  asymptotic  increase  in 
roll  moment  and  stability.  In  the  physical  sense,  however, 
the  vortex  may  not  be  capable  of  moving  so  close  to  the  fin 
as  to  cause  the  roll  moment  to  approach  infinite  values. 

The  interaction  between  the  vortex  and  its  image  results  in 
the  vortex  (and  image)  moving  outboard  and  parallel  to  the 
fin.  (The  reader  is  referred  to  Reference  6  for  a  theoreti¬ 
cal  treatment  of  this  interaction.)  Also,  the  sign  of  the 
rolling  moment  in  the  vicinity  of  a  fin  depends  on  which 
side  of  the  fin  the  vortex  is  located.  In  Figure  8,  for 
example,  the  sign  convention  of  the  roll  moment  indicates 
that  the  vortex  attracts  the  fin  as  it  approaches  it. 

Figure  9  shows  the  sensitivity  of  the  roll  moment 
to  perturbations  made  in  vortex  position  (I2.1A  ) .  This 
sensitivity  is  most  pronounced  when  vortex  1  lies  in  close 
proximity  to  the  vertical  fin.  In  the  vicinity  of  the  ver¬ 
tical  fin  tip,  the  roll  moment  becomes  highly  stable  when 
vortex  1  is  to  the  left  of  the  fin  and  highly  unstable  for 
positions  to  the  right  of  the  fin.  (Refer  to  Figure  7.) 

The  influence  is  less  pronounced  at  roll  angles  where  the 
vortex  is  located  further  away  from  the  fin.  A  similar 
situation  exists  in  Figure  10  where  perturbations  in  posi¬ 
tion  are  made  for  vortex  2.  Here,  the  roll  moment  shows 
greatest  sensitivity  when  vortex  2  is  in  close  proximity 
to  the  horizontal  fin.  Included  in  this  figure  is  the 
effect  of  simultaneously  increasing  the  perturbations  of 
both  vortices  at  =  17.9°. 

The  variation  of  roll  moment  with  fin  span  is 
shown  in  Figure  11.  It  is  apparent  that  large  moments  exist 


whenever  a  body  vortex  is  located  close  to  a  fin  tip.  In¬ 
creasing  the  span  reduces  the  roll  moment  to  an  asymptotic 
value  which  indicates  that  the  unbalanced  loads  inducing 
the  roll  are  located  near  the  inboard  portions  of  the  fins. 
It  follows  that  for  any  particular  condition  of  body  vor¬ 
tex-fin  geometry,  a  particular  value  of  span  will  result 
in  zero-induced  rolling  moment. 

An  almost  linear  variation  of  roll  moment  with 
perturbations  in  vortex  strength  is  indicated  in  Figures 
12  and  13.  Again,  the  greatest  variation  in  roll  moment 
occurs  for  a  vortex  position  close  to  the  fin. 

Results  of  Figures  8-13  are  cross-plotted  in 
Figures  14-20.  These  figures  are  felt  to  be  self  explana¬ 
tory  based  on  previous  discussions  and  are  therefore  not 
described  in  detail. 

Before  analytical  work  is  continued  by  including 
unsteady,  viscous  and  compressible  effects,  it  is  strongly 
suggested  that  measured  wind  tunnel  data  be  obtained  to 
support  the  results  of  this  report.  Further,  since  Refer¬ 
ence  3  does  not  consider  stagnation  point  locations  on  the 
body,  it  is  suggested  that  actual  locations  of  the  stagna¬ 
tion  point  be  determined  from  experiment  and  compared  with 
the  theoretical  locations  assumed  in  Appendix  B. 
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APPE1IDIX  A 


Derivatives  of 

The  derivatives  of  ^C'5)  are  listed  here  for  reference. 
From  Equ .  (5 } 

-/(=)=  (iVc^C^-c’)  ’ 

where 


Hence 


^  +  7i^ 


:W 


where 


_*»_  <^)|i 

'  (?vc,0-¥^s^+(^+c^2(rvc^vc%'f+c*ii-‘f/?'s» ' ' 


Also 


3^±C1  c,'(t) 


(?*+£:*)(%*-  c’) 
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APPENDIX  B 


The  purpose  of  this  appendix  is  to  find  an  analytical 
expression  for  the  free  vortex,  ¥  ,  as  well  as  the  location 
of  the  stagnation  point  to  satisfy  the  Kutta  condition  for 
the  flow  described  by  Equ.  (3).  The  limited  case  will  be 
considered  where  P  P2=0* 

To  find  the  following  procedure  is  used:  The 
stagnation  point  is  assumed  to  "move"  with  the  cross  flow 
vector,  0,  as  it  is  permitted  to  rotate  about  the  body. 
Thus,  in  the  Z  plane  the  stagnation  point,  Zq,  is  located 
on  the  body  contour,  as  follows, 

I  -tb  cf>=0 

Z0=  {  Zel(jr+0'>  o«P<tt/2 

[-ib  (p-Tr/T. 


To  obtain  the  corresponding  point  in  the  plane,  we 
make  use  of  the  conformal  mapping  (see  Equ.  (2a)) 


Substitution  of  (59)  into  (60)  gives,  for  0  «$><7T/Z 
y  e-Y R*COS  2 (<t> i-ir)  ±1  Vc*- Z«Cos*Z(&+Tr)  . 


(60) 


Tie  complex  velocity  is  obviously  zero  at  a  stagnation 
point.  Therefore,  if  we  first  differentiate  Equ.  (3)  to 
obtain  (.the  complex  velocity)  and  equate  this  result 

to  zero,  we  have 
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After  algebraic  reduction  we  obtain  the  general  ex¬ 
pression 

- 

Y &£')’*  2MT  &+ -ff )  sin  0 +  -§r)cos0 

«■ 

-r,(i  + 

+r2(i+ 

Obviously,  for  the  simple  limited  case  where 

fj-  =t  ~  O  ,  we  have 

r- 

Y(<t>,f£)=2Trxr  (?'+•§!  )smip  +  l^!r^r)cos<p. 

%  - 

To  express  Y  as  a  real  quantity,  we  can  further  define 


Substitution  in  the  above  expression  gives 

L}-rrTJc(coso0sir\<p  -sine0coscp^  m 
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It  is  easy  to  show,  by  substitution*  that  V  is  zero  wher¬ 
ever  flow  symmetry  exists  about  the  configuration;  that  is 
where  0,  77/ 4  and  77/2- 


